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ABSTRACT 

We present the first application of the 3D cosmic shear method developed in Heavens 
et al. (2006) and the geometric shear-ratio analysis developed in Taylor et al. (2006), 
to the COMBO- f 7 data set. 3D cosmic shear has been used to analyse galaxies with 
redshift estimates from two random COMBO-17 fields covering 0.52 square degrees in 
total, providing a conditional constraint in the (as, Q m ) plane as well as a conditional 
constraint on the equation of state of dark energy, parameterised by a constant w = 
Pde/pdeC 2 . The (os, fi m ) plane analysis constrained the relation between as and fl m to 
be crg(ri m /0.3)°' 57±019 = 1.06lg *g, in agreement with a 2D cosmic shear analysis of 
COMBO-17. The 3D cosmic shear conditional constraint on w using the two random 
fields is w = — 1.27^0 70- The geometric shear-ratio analysis has been applied to the 
A901/2 field, which contains three small galaxy clusters. Combining the analysis from 
the A901/2 field, using the geometric shear-ratio analysis, and the two random fields, 
using 3D cosmic shear, w is conditionally constrained to to = — 1.08lgg|. The errors 
presented in this paper are shown to agree with Fisher matrix predictions made in 
Heavens et al. (2006) and Taylor et al. (2006). When these methods are applied to 
large datasets, as expected soon from surveys such as Pan-STARRS and VST-KIDS, 
the dark energy equation of state could be constrained to an unprecedented degree of 
accuracy. 
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1 INTRODUCTION 

This paper presents the first application of 3D weak lens- 
ing techniques developed in Heavens (2003), Heavens et al. 
(2006), Jain & Taylor (2003) and Taylor et al. (2006) to data. 
The data used is the COMBO-17 survey (Wolf et al. 2001; 
Wolf et al., 2004) which is a multi-band photometric survey 
with exceptional image quality and is ideal for a 3D weak 
lensing study. The power of these methods is cosmological 
parameter estimation, focussing especially on measuring the 
equation of state of dark energy which appears to be respon- 
sible for the acceleration of the Universe. 

The case that the isotropic expansion of the Universe is 
accelerating is now convincing. The acceleration can be at- 
tributed to the effect of a negative pressure component, dark 
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energy, which accounts for approximately 70% of the mass- 
energy of the Universe. However the identity of dark energy 
is entirely unknown. The nature of dark energy may be ul- 
timately determined by establishing its equation of state, 
parameterised by 

W = Pdc/pdcC 2 . (1) 

A cosmological constant has w = — 1, a dynamical dark en- 
ergy such as quintessence may have w / -1. Note that in 
general w may be a function of redshift, z. This study at- 
tempts to use 3D weak lensing to constrain w (where we 
have assumed that w is a constant), for the first time. It 
is based on only 0.78 square degrees of data, and is essen- 
tially a proof of concept in preparation for much larger sur- 
veys such as VST-KIDS, Pan-STARRS (Kaiser, 2005) or the 
Dark Energy Survey (Wester, 2005) which could lead to very 
accurate measurements of w and its redshift evolution. 

We also present constraints on the amount of matter 
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in the Universe £l m and the clustering of matter, parame- 
terised by as, the rms of the fractional mass density fluc- 
tuations in spheres of radius 8/i -1 Mpc. Weak lensing has 
already proven to be a powerful probe of both fi m and as, 
using 2D weak lensing techniques. We show that a fully 3D 
shear analysis can also place tight constraints on the matter 
content and clustering. 

Other surveys have used weak lensing data to con- 
strain cosmological parameters using 2D and tomographic 
tests. Most recently Semboloni et al. (2006), Hoekstra et 
al. (2006), Schrabback et al. (2006) and Hetterscheidt et al. 
(2006) have all constrained as, 0. m and fl de , though all these 
surveys cover a much larger area than COMBO-17 . Hoek- 
stra et al. (2006) use the first data release of the CFHTLS 
Wide survey, which covers 22 square degrees in the i'-band, 
to constrain a 8 = 0.85±0.06 for ft m = 0.3, the full CFHTLS 
survey will cover 170 square degrees in 5 photometric bands. 
Semboloni et al. (2006) use the CFHTLS Deep data which 
covers 2.34 square degrees in (u* , g' , r' , i' , z') to constrain 
as = 0.89 ± 0.06 for Q m = 0.3. Semboloni et al. (2006) also 
combine with Hoekstra et al. (2006) to place an upper bound 
on w, marginalizing over O m , of w < —0.8. 

COMBO-17 has the best and most reliable photometric 
redshifts to date, due to the large number of bands, and so 
is ideal as a survey to test the 3D weak lensing constraints. 
COMBO-17 is then ideal for this proof of concept, however 
as shown in Heavens et al. (2006) and Taylor et al. (2006) 
when much larger survey areas are available, a 5-band large 
area survey could constrain the dark energy equation of state 
much better than correspondingly smaller area 17-band sur- 
vey. 

The first 3D weak lensing technique used in this pa- 
per is a 3D cosmic shear analysis which analyses the full 
3D cosmic shear field using a spherical harmonic expansion, 
proposed in Heavens (2003) and developed in Castro et al. 
(2005) and Heavens et al. (2006). The second technique, sug- 
gested by Jain & Taylor (2003), is a geometric shear-ratio 
analysis which takes the ratio of tangential shears around 
galaxy clusters, developed in Taylor et al. (2006). Both of 
these analyses are applied to the data available. We use 3D 
cosmic shear to place conditional constraints on w and to 
place conditional constraints in the (as, O m ) plane using the 
two random fields. The A901/2 field is centered on a known 
supercluster, so we use the geometric test on this field to 
place a conditional constraint on the dark energy equation 
of state parameter w. This separation of the data allows for 
a combination of the two methods in the constraint of to. 

The results of this paper are a proof of method for these 
techniques. Heavens et al. (2006) and Taylor et al. (2006) 
show that in order to constrain the dark energy equation of 
state to Aw « 0.01 large and deep photometric surveys will 
be needed. The errors on the results in this paper are com- 
pared to predictions made using the Fisher matrix formal- 
ism used in Heavens et al. (2006) and Taylor et al. (2006). 
Brown et al. (2003) have already applied 2D weak lensing 
to the COMBO-17 data set to constrain the (as, f2 m ) plane. 
The results presented in this paper using the fully 3D cosmic 
shear analysis on the same data set, should be in approxi- 
mate agreement (but do slightly better than) a 2D analysis. 
It should also be noted that in order to constrain the dark 
energy equation of state one necessarily needs 3D methods. 
To determine whether dark energy is a field or a manifesta- 



tion of of modified gravity (e.g. Ishak et al., 2006), one needs 
methods can probe the expansion history of Universe, as is 
the case with the geometric shear-ratio analysis, or both the 
expansion history and the growth of structure as is the case 
in the 3D cosmic shear analysis. 

The structure of this paper is as follows: firstly the 
COMBO-17 data set will be introduced and discussed in 
Section 2, the application of the 3D cosmic shear analysis 
to the two random fields will be presented in Section 3, the 
geometric shear-ratio analysis and the results of applying 
the method to the A901/2 field will be presented in Section 
4. The constraints from the A901/2 field (using the geomet- 
ric shear-ratio analysis) and the constraints from the two 
random fields (using the 3D cosmic shear analysis) will be 
combined in Section 5. Conclusions will be presented in Sec- 
tion 6. 



2 THE COMBO-17 SURVEY 

The COMBO-17 survey is a 17-band photometric redshift 
survey with gravitational lensing quality R-band data (Wolf 
et al., 2001; Wolf et al., 2004). The survey consists of five 
fields each covering 0.26 square degrees. All of the fields 
were observed using the Wide-Field Imager (WFI) at the 
MPG/ESO 2.2m telescope on La Silla in Chile, with a 4 x 2 
array of 2048 x 4096 pixel CCDs, each pixel subtending 0.238 
arcseconds. 

In this paper we will use three of the COMBO-17 fields, 
which were observed and reduced earlier than the remaining 
two, and for which there are redshift estimates and a shear 
catalogues available. One of the fields used, the A901/2 field, 
is centred on the Abell 901/2 supercluster which has previ- 
ously been analysed in 2D by Gray et al (2002) and in 3D 
by Taylor et al. (2004). The A901/2 supercluster consists of 
three smaller clusters; A901a, A901b, A902, all at a redshift 
of z w 0.16. It should be noted that supercluster refers to 
a 'web of clusters', the individual clusters are much smaller 
~ 10 14 M Q (see Taylor et al., 2004) than large strong lens- 
ing clusters for example A 1689. For an individual cluster the 
fractional error on w should decrease as the mass of the clus- 
ter increases. However in a large area survey there should 
be many more low and medium mass clusters than large 
clusters so that the constraint on w is dominated by the nu- 
merous medium mass clusters; for a detailed discussion see 
Taylor et al. (2006). 

The COMBO-17 survey also observed a randomly se- 
lected area of sky, and a relatively empty, but well observed 
area. The Sll field was a randomly selected area of sky, 
that contains a moderately large cluster Abell 1364 at a 
redshift of z ~ 0.11. The CDFS field was chosen to overlap 
the Chandra Deep Field South, a relatively 'empty' region of 
sky containing no significant galaxy clusters. We only used 
galaxies with reliable photometric redshifts and with an R 
magnitude of R < 24. 

2.1 Photometric Redshifts 

Each of the COMBO-17 fields was observed in 17 different 
filters, with the intention of obtaining object classification 
and accurate photometric redshifts. In order to provide reli- 
able redshifts, the filter set included five broad-band filters 
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(UBVRI) and 12 medium-band filters from 350 to 930 nm. 
This observing strategy allows simultaneous estimates of 
Spectral Energy Distribution (SED) classifications and pho- 
tometric redshifts from empirically-based templates. Wolf et 
al. (2001) describe in detail the photometric redshift estima- 
tion methods used to obtain typical accuracies of a z « 0.05 
for galaxies throughout < z < 1. We fit an empirical line 
to the data such that 

o z (z) =0.03(1 + 2) 1 ' 5 , (2) 

this <r z (z) is used in the likelihood analysis. It should be 
noted that the parameter constraints are not sensitive to 
the exact functional form of a z (z). 



2.2 Shear Measurements 

Throughout the observing campaign the R filter was used 
in best seeing conditions, in order to provide a deep 7?-band 
image from which to measure the gravitational shear. Gray 
et al. (2002) discuss the procedure used to reduce the R band 
imaging data, which totalled 21 hours for the three fields 
used. As described by Gray et al. (2002) and Brown et al. 
(2003) the 352 individual chip exposures for each field were 
registered using linear astrometric fits, with a 3a rejection 
of bad pixels and columns. 

The Kaiser, Squires & Broadhurst (1995; KSB) weak 
lensing measurement method was applied, using the imcat 
shear analysis package, to the reduced images (see Gray et 
al., 2002; Brown et al., 2003). This resulted in a catalogue 
of galaxies with centroids and shear estimates throughout 
the fields, corrected for the effects of anisotropic smearing 
and point spread function (PSF) circularisation. The pho- 
tometric redshift estimates were appended to this catalogue 
for each galaxy from the standard COMBO- 17 analysis of 
the full multi-colour dataset. Of the 37,243 galaxies in the 
shear catalogue, 36% have a reliable photometric redshift, 
the remainder being fainter than the R = 24 reliability limit 
of the redshift survey. The requirement for the 3D lensing 
study, that the redshift of each galaxy be known, clearly re- 
sults in an immediate reduction of available galaxies. It is 
apparent that most of the background sample is composed 
of galaxies that are small, and fainter than the magnitude 
limit of the redshift survey. These catalogues are the raw 
data used in this analysis. Brown et al. (2003) also include 
galaxies without assigned redshifts into their analysis, this 
can also be done in the case of 3D weak lensing however 
since this paper is a proof of concept for the 3D weak lens- 
ing methods the galaxies without redshifts will be left out 
of this analysis. 



3 THE 3D COSMIC SHEAR ANALYSIS 

We have applied the 3D cosmic shear method to the CDFS 
and Sll fields of COMBO-17 in order to constrain w and 
jointly constrain as and Q m . Bacon et al. (2005) analysed 
COMBO-17 using a real-space 3D cosmic shear method to 
constrain the evolution of dark matter clustering. The re- 
sults presented in this Section are based on the methods 
outlined in Heavens et al. (2006). 



3.1 3D Cosmic Shear Likelihood 

The statistics we choose to work with are the transforms of 
the galaxy ellipticities. The transforms are are defined, for a 
given radial fc-mode and angular £-mode, by summing over 
all galaxies g, each at a redshift z and angular position g , 
in a given field catalogue 

7i(M) = y|^effcj £ (fc^)^e-^^M/( z ). (3) 
s 

The hat indicates that these are estimators of the transform 
of the shear field, but the relationship is not direct as the 
number density of sources is non- uniform. 

g refers to the galaxies in the sample, ef is the i com- 
ponent of the complex ellipticity of the galaxy related to 
the weak shear 7? and the intrinsic ellipticity of the galaxy 
by e i = e i (intrinsic) + 7?. r° g denotes the comoving dis- 
tance to a galaxy calculated from the photometric redshift 
of the galaxy by assuming a fiducial cosmology. The ji(z) 
are spherical Bessel functions. W(z) is a weighting function 
which we set to W(z) = 1 for the remainder of this paper; 
for an investigation of the effect of changing the weighting 
scheme see Heavens et al. (2006). The Xg factor is intro- 
duced as a result of relating the potential to the shear field 
j(r) = i S (5<j>(r) ( see Heavens et al., 2006) and is given by 

v — (^y ~ ^x) + "Zif-xZy ... 
X £ = J2 • ( 4 ) 

The expansion using spherical Bessel functions is natu- 
ral for a flat universe. For non-flat universes, the appropriate 
functions are ultra-spherical Bessel functions $^(t/), but in 
the i 3> 1 and k S> (curvature scale) _1 regime these are well 
approximated by ordinary Bessel functions &p(y) — > je(kr) 
(Abbott and Schaefer, 1986; Zaladarriaga and Seljak, 2000). 

For each ef component there is a real and imaginary 
estimator (via the e~^ 8g factor) so that the whole data 
vector used in the likelihood analysis consists of four inde- 
pendent vectors at each k and £ : 7P, 7*, 7^, 7I, where the 
R superscript denotes the real part of the 7, estimator and 
I the imaginary part. Note that this is for a given £-mode. 

The fiducial cosmology, denoted by the superscript 0, is 
in this case chosen to have Q m = 0.3, fi de = 0.7, Qb = 0.04, 
h = Ho/100 kms^Mpc -1 = 0.71, ct 8 = 0.8, w = -1.0; 
we also set the scalar spectral index to be n s = 1.0 and 
its running a„ = 0.0. The choice of the fiducial cosmol- 
ogy does not affect the results presented in this paper, as 
long as the measured shear estimates and theoretical cal- 
culations use the same fiducial cosmology to calculate the 
transform coefficients. This fiducial cosmology simply acts, 
via the spherical Bessel functions, to weight the shear val- 
ues in a particular way. The cosmological dependence comes 
from the shear values themselves, 7?, the cosmological de- 
pendence of the calculated covariance matrices come from 
modelling the shear-shear covariance. We tested a variety of 
fiducial models and the results were indeed unaffected. 

Note that the average value of -yi(k,£) is zero, so that 
information on the cosmological parameters comes from the 
dependence of the signal part of the covariance matrix C 
i.e. we adjust the parameters until the covariance of the 
model matches that of the data. This was the approach of 
Heavens and Taylor (1995); Ballinger, Heavens and Taylor 
(1995); Tadros et al. (1995); Percival et al. (2004) in analy- 
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sis of large-scale galaxy data. The details of the covariance 
matrix derivation are given in Heavens et al. (2006), where 
the covariance matrix is given as the sum of signal and noise 
terms C = S+N. The signal part of the covariance of y(k, £) 
for a survey of size AO x AO can be written as 

S = <7(M)7*(fc',0)s 



= I ^Q(e,£',lk,k')\X- £ \' 

r A6/2 . ~. r Aff/2 

/ <fOe-« 1 -*- / 
J-ab/2 J-Ae/2 



The Q matrix is given by 



d 2 , e -i ( t'-£).e< 



(5) 



Q (e,e',lk,k') = ^0- J ^G(£,lk,k)G(£',l,k',~k) 



(6) 



~£,k,k)=k J dz dz p n z (zp)p(z p \z)U(£,r,k)ji(kr°) 



where 
G{£ 

(7) 

where z v is an integral over redshift given an assumed cos- 
mology, and the integral over z uses the cosmology to be 
tested. n z (z) is the predicted number density of objects as 
a function of redshift which is measured from the survey 
in question. p(z p \z) is a probability distribution in redshift 
which, by convolving with the redshift distribution, takes 
into account the uncertainty in redshift. We assume a Gaus- 
sian probability distribution with the width being the mea- 
sured photometric redshift as a function of redshift (see Sec- 
tion 2.1). 

The U matrix used in equation (7) is 
F K (r, f) 



U(£,r, k) 



f 

Jo 



df ■ 



a(f) 



-y/P^Fjjeikr), (8) 



where Ps(k;r) is the matter power spectrum for the cos- 
mology to be tested. We compute the nonlinear power spec- 
trum using the fitting formulae of Smith et al. (2003) , based 
on linear growth rates given by Linder and Jenkins (2003). 
Ffc(r,r') = (1/r' — 1/r) for a flat Universe (assumed in this 
paper) and a(r) is the dimensionless scale factor. 

The integrals over 6 in equation (5) can be evaluated 
so that 

S = J -^Q(£,£',£,k,k')\X £ \ 2 F(£,£',l), 

where 

T(£, £',£)= n - ^= — sin \{£-£)i^ 6 



(9) 



{£ - £)i{£ - £')i 



(10) 



i — x,y represents the x and y components of a Cartesian 
coordinate system for a survey. 

The shot noise part of the covariance matrix is calcu- 
lated by assuming a Poisson sampling an underlying smooth 
density field and is given by 

N = (%(k,£)%(k',£')) S N 

= ~t^T J dzn z (z)kj e (kr°)k'j e (k'r°)S^ l3 Su, (11) 

where Af2 = AO x A9. u e is measured from the data, for the 
CDFS field a € = 0.19 and for the Sll field <r £ = 0.22. 



For large surveys the calculation simplifies as T — > 
(A9) 2 8fg8f,g (see equation 10) and the covariance matrix 
becomes diagonal in (£, £'). For COMBO-17, the survey is 
small so that it is necessary to compute the integral in equa- 
tion (9) accurately using the full J-{£, £',£). We do this for 
the diagonal components of Q, but we make an approxima- 
tion by ignoring the off-diagonal components. Improvement 
on this would involve computing the full covariance resulting 
in vast computational expense, which is not really warranted 
by the size of the dataset. All correlations between fc-modes, 
for any given £-mode are fully taken into account. 

We assume then that the distribution of ^-modes can 
be represented by a multivariate Gaussian. The likelihood 
function for a given £ and set of cosmological parameters 
{0 a } is given by 



-2\n L e (0 a \D) 



A={R,I} i={l,2} 

+ H\C i Ak,k') AA \) 

+J2rtmc- e l (k,k')) 

kk> 



(12) 



AA t' T (k') 



where A = {R, 1} is a sum over the real and imaginary 
estimators and i = {1,2} is a sum over the 71 and 72 shear 
components. N A is the number of fc-modes in the N A x N A 
covariance matrix Ci/(fe, k') AA . The log-likelihood is then 
summed over each independent £ — (£ x ,£y) mode. Note that 
since we have four independent data vectors, two real and 
imaginary pairs, we only investigate the range l x > to 
avoid double counting. 

The I and k ranges and resolutions used are as fol- 
lows. In integrating over £ in equation (5) we found that 
the signal converges at Ali — 100 and for the range 
(£i - 1500) < h < [£i + 1500) where i = x,y. The k res- 
olution used was Afc = 2 x 10~ 3 Mpc -1 , as it was found 
that the signal part of the covariance matrix converges at 
Afc w (2-7r/r m ax) ~ 2 x 10 -3 Mpc -1 where r max is the dis- 
tance corresponding to a maximum redshift of 2 ~ 1. The fc 
range used was 0.01 < k < 1.5 Mpc -1 . The I values avail- 
able are constrained by the survey geometry, ii = 2gi where 
n is an integer; i\ = 671. We tested the lower £ limit and 
found no change in the cosmological constraints by using 
£1 — 700 instead, showing that these results are robust to 
the details of the lower i range used. We use all modes 
with \£\ < 2500 to avoid the highly non-linear regime in 
which baryonic physics may have a significant effect (Zhan 
& Knox, 2004; White, 2004) on the power spectrum. We cal- 
culate non-linear effects accurately using the fitting formula 
of Smith et al. (2003) for the matter power spectrum. The 
use of a Gaussian likelihood is an approximation. The scales 
probed by the technique currently extend significantly into 
the non-linear regime, but the integral nature of lensing is 
expected to make the shear more Gaussian (Hu & White, 
2001). However, it remains to be tested with simulations 
whether the Gaussian approximation leads to significant er- 
ror. This will be the subject of future work. Note that since 
£ = (£ x ,£ y ) the £ range used corresponds to 26 independent 
£- modes. 

The use of spherical Bessel functions in the coefficients 
used means that, for any given £-mode, there is a range of fc 
for which the signal is zero up until a particular value of fc « 
|^|/r ma x (see Castro et al., 2005). These zero modes result in 
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Figure 1. The solid lines show the two-parameter 1-<t conditional 
constraints in the (as, f2 m ) plane from applying the 3D cosmic 
shear analysis to the CDFS and Sll fields only. The dashed con- 
tours show the two-parameter l-cr conditional constraints from 
the Brown et al. (2003) analysis using 50% more fields: CDFS, 
Sll and A901/2. 



singular covariance matrices, however this behaviour can be 
taken into account using the prescription given in Appendix 
A. 



3.2 3D Cosmic Shear Results 

This Section presents the result of applying the 3D cos- 
mic shear analysis to the CDFS and Sll fields. The results 
will be compared with the Fisher matrix analysis of Heav- 
ens et al. (2006) and with the 2D cosmic shear analysis of 
Brown et al. (2003). Unless otherwise stated the fiducial 
cosmology that will be assumed throughout this Section is 
n m = 0.3, n de = 0.7, n h = 0.04, h = 0.71, a 8 = 0.8, w = 
— 1.0, n s = 1.0, a„ — 0.0, any constraints for particular pa- 
rameters are conditional on these values. 

Figure 1 shows the two-parameter l-cr contours from 
applying the 3D cosmic shear analysis to the CDFS and 
Sll fields only. The dashed line in Figure 1 shows the two- 
parameter 1-ct contours from Brown et al. (2003) where a 
traditional 2D cosmic shear analysis was performed on all 
three COMBO-17 fields, CDFS, Sll and A901/2 using only 
galaxies with accurate redshifts. It can be seen that the 3D 
cosmic shear analysis constrains a very similar area in the 
(as, Qm) plane, particularly at the concordance values of as 
and Q. m using less than two thirds the number of galaxies 
used in the 2D analysis (since the A901/2 field contains more 
galaxies than the CDFS and Sll only 63% of the galaxies 
used in the 2D analysis have been analysed). 

A common way to express the constraint in the (as, ii m ) 
plane is to constrain the parameterisation ag^m/OS) 13 = Q, 
where f3 expresses the curvature of the constraint and a the 
normalisation of the curve. The 3D cosmic shear analysis 



constrains these parameters to be 
a = 1.06t° .\l 

(3 = 0.57±C;11 (13) 

these are consistent with the constraints from the 2D anal- 
ysis of Brown et al. (2003). It should be noted however that 
Brown et al. (2003) have a maximum I limit of ~ 10, 000, so 
that they use substantially more modes in the angular direc- 
tion. Also, their main result of CT 8 (^m/0.3) ' 49 = 0.72±g;g| 
included galaxies with unknown redshifts, the result shown 
(dashed line in Figure 1) is their result when considering 
galaxies with only reliable redshift estimates. Since the range 
in £-modes and the number of galaxies are different in the 
two analyses this comparison cannot make any conclusions 
on the relative merit of the two techniques. It is sufficient to 
say that our results agree with Brown et al. (2003) when we 
include the same sort of data i.e. from an analysis in which 
the majority of the galaxies were the same the results are 
consistent with one another. 

The high value of erg is unexpected for the CDFS and 
Sll fields, the CDFS result on it own favours a lower as 
(with substantially increased error contours), so that the 
cluster in the Sll field appears to increase the most likely 
value of as- Also the effect of not including galaxies with 
photometric redshifts appears to increase the most likely 
clustering value in Brown et al. (2003) (compare Figures 19 
and 20 in Brown et al., 2003). Since we do not include galax- 
ies for which a redshift is unknown we may expect to find a 
high clustering value in a similar way to Brown et al. (2003) . 
We only used galaxies with reliable redshifts since this paper 
is a proof of concept for the 3D weak lensing methods, how- 
ever galaxies with unknown redshifts could also be included 
in this analysis. 

The Fisher matrix calculations in Heavens et al. (2006) 
can be used to predict the estimated uncertainties from 
this analysis. Fisher matrix predictions, by construction, 
predict Gaussian likelihood surfaces, the curved constraint 
shown here highlights one limitation of the Fisher matrix 
technique to predict uncertainties when the errors are so 
large. However, using the techniques outlined in Heavens 
et al. (2006), we predict, for two COMBO-17 fields, a con- 
ditional constraint of Aas = 0.19 (assuming Q. m = 0.3). 
This is in agreement with the measured conditional error of 
as(Qm = 0.3) = 1.05±0.20, highlighting that the predictions 
made in Heavens et al. (2006) are reliable. The values used 
in the Fisher matrix calculation were an area of A = 0.52 
square degrees, to a median redshift of z mc dian = 0.8 using 
the photometric redshift error given in equation (2). 

Figure 2 shows the conditional constraint on w from the 
CDFS and Sll field only using the 3D cosmic shear analy- 
sis. The constraint is asymmetric in that the range w < — 1 
is more likely than w > — 1. This is due to the fact that 
values of w < —1 represent dark energy scenarios in which 
the dark energy density is less in the past, so it is more 
difficult to constrain its equation of state. Semboloni et al. 
(2006) also found a similar asymmetric constraint when us- 
ing weak lensing tomography applied to the CFHTLS sur- 
vey. The conditional constraint on w is 

™ = -1.27lE , ) ; 6 70- (14) 

This result is consistent with other observations (for example 
Spergel et al., 2006) and with a cosmological constant model 
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Figure 2. The one-parameter maximum likelihood constraint on 
w from the CDFS and Sll fields using the 3D cosmic shear anal- 
ysis. The dashed line shows the most likely value and the dot- 
dashed show the one-parameter l-cr constraints. 



for dark energy. The Fisher matrix calculations, presented 
in Heavens et al. (2006) predict a conditional error on w 
from two COMBO-17 fields to be Am = 0.62 which is in 
agreement with the constraints presented here. 

Typical reduced \ 2 values for a given ^-mode in the 
CDFS and Sll fields analyses are Xcdfs ~ 1-01 and Xsu ~ 
0.98, the number of degrees of freedom for a given £-mode 
are the corresponding number of non-singular fc-modes used 
in the analysis, typically ~ 600 for an average £-mode. The 
range of \ 2 values are consistent with a good fit to the data. 



4 THE GEOMETRIC SHEAR-RATIO 
ANALYSIS 

We have applied the geometric shear-ratio analysis to the 
A901/2 field of the COMBO-17 survey in order to condi- 
tionally constrain w, and compare the measured constraint 
with the predicted constraint from a Fisher matrix calcula- 
tion. The results presented in this Section are an extension 
and use of the methods outlined in Taylor et al. (2006). 



4.1 Geometric Shear-Ratio Likelihood 

To implement the geometric shear-ratio analysis, we first se- 
lected the peaks in the convergence field of the three clusters, 
A901a, A901b and A902. Taylor et al. (2004) have shown 
that there is a fourth cluster, CB1, in this field, which lies 
behind A902 at a redshift of z = 0.42. Here we shall ignore 
the contribution of this cluster, although this will in princi- 
ple bias our results slightly. To estimate the effect of the bias 
the CB1 cluster increases the tangential shear, at z > 0.4, by 
5~ft < 0.02 (see Taylor et al., 2004). Using the simple error 
formula from Taylor et al. (2006) this increase in tangential 
shear may bias the value of w by Sw < +0.03. 

We use the positions of the centre for each cluster given 
by Taylor et al. (2004) , in which the tangential shear around 



each cluster is used to determine the 3D position and mass 
of the clusters. We averaged the tangential shear in annuli 
around each cluster in a series of redshift bins, following 
Taylor et al. (2004; see Figure 3), the width of the redshift 
bins is equal to the photometric redshift error at the redshift 
of the bin using the result from Section 2.1. The lensing 
signal from a cluster is given by the tangential shear 



7t 



-[71 cos(2(^) + 72 sin(2^)], 



(15) 



for a given redshift bin. The error on the tangential shear 
was estimated by the orthogonal shear signal 



7x = [-71 sin(2ip) + 72 cos(2y)]. 



(16) 



Since a lensing cluster induces a tangentially aligned shear 
signal, any orthogonal component is assumed to be due to 
noise. In this case we define a polar coordinate system about 
the cluster centre and assume a circularly-symmetric mass 
profile, which should produce a purely tangential signal. Un- 
der these assumptions any orthogonal component is taken 
to be due to measurement error, as opposed to deviations 
from circular symmetry in the lens (which can also produce 
an orthogonal component in the circular coordinate system 
chosen about the cluster center). The tangential shear in 
each angular and redshift bin was then fitted with a least- 
square fit to a singular isothermal sphere (SIS) profile; 



7*,sis(, 



(17) 



8e{z) is the Einstein ring radius which parameterises the 
amplitude of the tangential shear as a function of source 
redshift. Note that the assumption of a SIS is not necessary, 
in the case of a large data set the average tangential shear in 
an aperture could be measured directly, and the orthogonal 
shear component for the error, with no assumption on the 
radial tangential shear profile made. Since this data set con- 
sists of only three small clusters the SIS was adopted so that 
a signal could be measured, and for the radii from the centre 
of the clusters probed should be an adequate approximation. 

Here D denotes data and R is the theoretical estimate 
for the shear ratio, dependent on cosmology. The theoretical 
ratio of shears, Rij for a pair of redshift bins could then be 
estimated by 



R _ 0E(zj) _ (riens -r i )/r l 
13 E (zj) (nens - rj)/rj 



(18) 



where r is the predicted comoving distance for a given cos- 
mology. The data is simply the ratio of tangential shears 



lt(Zi) 



Dij = 



(19) 



The measured ratio Dij and the calculated ratio Rij are then 
used in the likelihood function, summing over all pair-pair 
configurations, given by equation 

-2lnL c (Q de ,D. m ,w,w a \D) = 

J2(R, - D I1 )[C^]-\R V - D v ), (20) 

A 1 ! v 

for a given cluster. The notation for a pair of background 
bins has been compressed as /i — and v — (m, n), and 
all degenerate pair-pair combinations have been accounted 
for. The likelihood functions for multiple clusters are mul- 
tiplied. We can write the covariance matrix for shear ratios 
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w 

Figure 3. The dark energy geometric shear-ratio analysis ap- 
plied to the supercluster Abell A901/2. The dashed line marks 
the maximum likelihood value, the dot-dashed lines show the one- 
parameter 1-cr limits. Note that the x-axis scale has been extended 
relative to Figures 2 and 3 to encompass the confidence limits of 
this analysis. 

as 

C™ = (AR V AR^). (21) 

The full covariance matrix includes shot noise and cosmic 
shear terms. For a full description see Taylor et al. (2006). 

4.2 Geometric Shear-Ratio Results 

The results are shown as a ID likelihood plot in Figure 3. 
The result is conditional on Q m = 0.30 , Qde = 0.70 and 
w a — 0.0. The dashed line is our measured most likely value, 
the dot-dashed lines are the one-parameter 1-cr (68%) con- 
fidence limits. Our constraint on w is: 

w = -0.11±\H. (22) 

The constraint again shows an asymmetry between the w < 
— 1 and w > ~ 1 regions for the same reason given in the 3D 
cosmic shears constraint, that w < — 1 represents a lower 
dark energy density in the past. The minimum \ 2 value is 
Xmin = 122 which is consistent with the number of degrees 
of freedom in the experiment. Given that z max w 2.0 and 
Az ~ 0.05 and we analyse iV c i uster = 3 clusters the predicted 

Xmin = (-Zmax/A^iVclustor « 120 SO that Xrcducod = 1-01 and 

should be xLiucod = 1 ± 0.12. 

The result is consistent with other constraints on w, and 
the confidence limits allow for most dark energy models. It 
should be emphasised that this constraint comes from only 
three small clusters. 

The Fisher matrix calculations, in Taylor et al. (2006) 
predict a conditional constraint on w of Aw = 1.10 for 
COMBO-17 which was created by assuming only three clus- 
ters at z = 0.16 with M c i us t cr = 1O 14 M0. The predicted 
conditional constraint is approximately the same as the mea- 
sured constraint, thus verifying the Fisher matrix method- 
ology. 




-2.5 -2 -1.5 -1 -0.5 



w 

Figure 4. The one-parameter maximum likelihood constraint on 
w obtained by combining the geometric shear-ratio analysis con- 
straint from A901/2 field and the 3D cosmic shear analysis con- 
straint from the CDFS and Sll fields. The solid line shows the 
combined constraint, the dashed line shows the 3D cosmic shears 
constraint shown in Figure 2, the dot-dashed line shows the geo- 
metric shear-ratios constraint shown in Figure 3. 

5 A COMBINED CONSTRAINT ON w 

Since the A901/2 field was analysed separately from the 
CDFS and Sll fields the geometric shear-ratio analysis con- 
straint can be combined with the constraint from the 3D 
cosmic shear analysis. Figure 4 shows the result of adding 
the constraints shown in Figures 2 and 3 from the CDFS and 
Sll fields and the A901/2 field respectively. The resulting 
conditional constraint on w is 

w = -1.08lS;£. (23) 

This result demonstrates the value of combining the two 
techniques that can analyse distinct parts of the data. A 
region of particularly high density such as the A901/2 field 
would raise issues of sample bias in a cosmic shear/spectral 
approach, however the geometric shear-ratio analysis neces- 
sarily needs such areas. The effect of adding the geometric 
shear-ratio analyses constraint is the most likely value be- 
coming more positive, and a slight reduction in the error. 
The most likely value of w = —1.08 is in complete agree- 
ment with other observations (for example Semboloni et al, 
2006; Spergel et al., 2006) and is close to the value of w ex- 
pected if dark energy is a cosmological constant. The caveat 
on this conclusion is that it is conditional on the other cos- 
mological parameters being fixed, and the error on w is still 
fairly large. Finally it is surprising, given the size of the 
error, that the maximum likelihood estimate is so close to 
w = — 1. 



6 CONCLUSION 

In this paper we have applied the 3D weak lensing methods 
introduced in Heavens (2003) and Jain & Taylor (2003) and 
developed in Heavens et al. (2006) and Taylor et al. (2006) 
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to data for the first time. We used the COMBO-17 data set, 
a multi-band photometric survey, ideal for a weak lensing 
study. We used three of the fields, CDFS, Sll and A901/2 to 
place conditional constraints on the equation of state of dark 
energy, w and the (as, Sim) plane. The size of the dataset 
is small so we did not expect accurate constraints and this 
paper is essentially a proof of concept for the methods. To 
this end we compute conditional constraints and compare 
with Fisher matrix predictions. A full analysis of a larger 
data set should, of course, quote marginal errors. 

The first method used was the 3D cosmic shear analysis 
which uses spherical harmonics to describe the fully 3D shear 
field. Applying the 3D cosmic shear analysis to the CDFS 
and Sll field we conditionally constrained an area in the (as, 
il m ) plane described by cr 8 (f2 m /0.3)°' 57±0 - 19 = l-06±g;^. 
This is in agreement with Brown et al. (2003) in which a 2D 
cosmic shear analysis was performed on 50% more data. The 
application of 3D cosmic shear conditionally constrained the 
equation of state of dark energy to w = — 1.27jlo'?o- 

The second method used was the geometric shear-ratio 
analysis, which takes the ratio of the tangential shear around 
galaxy clusters at different redshifts. We applied this anal- 
ysis to the A901/2 field which contains three small clusters 
conditionally constraining w = — O.lSlJ;^- 

Combining the constraint on w from the geometric 
shear- ratio analyses application to the A901/2 field and the 
3D cosmic shear constraint from CDFS and Sll we condi- 
tionally constrain w = — 1.08±q;||. For discussions on the 
relative merit of the two methods and varying observing 
strategies see Heavens et al. (2006) and Taylor et al. (2006) . 
These papers also discuss the effects of systematics, num- 
ber of observing bands used and the effects of the assumed 
fiducial cosmology. 

The constraints presented here do not improve much 
on our cosmological understanding, they are however in 
agreement with the currently accepted concordance model of 
Spergel et al. (2006). The constraint on w from such a small 
data set is encouraging and, with the warning and caveat 
that it is a conditional error, it is a forteriori consistent 
with dark energy being a cosmological constant. 

In order for these results to become more complete we 
could marginalize over an increasingly large cosmological pa- 
rameter set. However this would rapidly result in a loss of 
any constraint, for such a small survey area, due to degen- 
eracies between the parameters as shown in Heavens et al. 
(2006) and Taylor et al.(2006). The result errors agree with 
the Fisher matrix predictions and are a very reliable proof 
of concept for these methods. 

The agreement of the results presented here with the 
Fisher matrix predictions using the methods presented in 
Heavens et al. (2006) and Taylor et al. (2006) are a vali- 
dation of the Fisher matrix framework and an encouraging 
sign that the predictions made in these papers are robust 
and accurate. COMBO-17 was an ideal survey upon which 
to test these 3D weak lensing methods, however with much 
larger area surveys with fewer observing bands, such as VST- 
KIDS, Pan-STARRS (Kaiser et al., 2005) or the Dark En- 
ergy Survey (Wester, 2005), these techniques could constrain 
the dark energy equation of state to approximately 1% at a 
particular redshift. 
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APPENDIX A: REMOVAL OF SINGULAR MODES 

In this Appendix we show how any singular (k, I) modes can be removed from the covariance matrices used in the 3D cosmic 
shear likelihood analysis. 

We begin with a square covariance matrix C which can be decomposed using a standard singular value decomposition 
(SVD) into 

C = UWV T , (24) 

where W is a diagonal matrix that contains the singular values. Note that U and V are eigenvector matrices of CC T and that 
U" 1 = U T and V -1 = V T . Our covariance matrices are symmetric so that U = V in this case. Now consider one of our data 
vectors y^(k,£) represented by x which can be transformed to a new data vector y via 

y = Bx (25) 

where B can be any, not necessarily square, transformation matrix. A new covariance can then be defined 

C'i] = iViVj) = {B lk x k Bjixi) = BikBjiC'ki, (26) 

which implies that 

C' = BCB T . (27) 

The choice of transformation, in this case, is motivated by decomposing C using a Cholesky decomposition which yields 
C = UWU T = LL T , where L = UW 1/2 . So we use B = W~ 1/2 U' 1 where W~ 1/2 = W 1/2 except that the elements of the 
inverse W matrix 1/wi have been replaced with zero if (l/w») >(threshold) where the threshold represents machine precision 
(see Press, 1993; Numerical Recipes). The matrix B now contains a band of values below which zeros remove any singular 
modes from either the data vector or the covariance matrix via equations (25) and (27). 

The transformation is performed using a fiducial cosmology (the choice of this does not affect the results) to yield a 
transformation matrix B which is then used throughout, y and C replace x and C in the likelihood analysis. 
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